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A Cosmic Battery 

loannis Contopoulos ^ and Demosthenes Kazanas ^ 


ABSTRACT 

We show that the Poynting-Robertson drag effect in an optically thin 
advection-dominated accretion flow aronnd active gravitating objects generates 
strong azimuthal electric currents which give rise to astrophysically signihcant 
magnetic fields. Although the mechanism is most effective in accreting 
compact objects, it seems very promising to also account for the generation of 
stellar dipolar fields during the late protostellar collapse phase, when the star 
approaches the main sequence. 

Subject headings: MHD — Sun: magnetic fields 


1. Introduction 

The origin of cosmic magnetic helds remains one of the open questions in astrophysics. 
Their origin has been sought in the entire history of the Universe going as far back as the 
QCD phase transition (Sigl, Olinto & Jedamzik 1997) the electroweak phase transition 
(e.g. Baym, Bodeker & McLerran 1996; Joyce & Shaposhnikov 1997) or even close to the 
Planck time (e.g. Ratra 1992). Most theories attempt to produce the observed magnetic 
fields by starting with some ‘seed magnetic field’ and follow its amplification by means of a 
‘dynamo mechanism’, i.e. a mechanism which can convert kinetic energy of the conducting 
fluid, into which the magnetic field is supposedly frozen, into magnetic field energy density 
(Moffat 1978; Parker 1979). Unfortunately, when one considers the back reaction of the 
stretched and folded magnetic field on the dynamics of the conducting fluid, dynamo 
action seems to be very ineffective, in that as the held at the small diffusive scales reaches 
equipartition, the large-scale component remains several orders of magnitude weaker than 
astrophysical observed magnetic helds (Vainshtein & Cattaneo 1992). 

In addition to these problems with the ehectiveness of the dynamo theory, the origin of 
the ‘seed magnetic held’ needed to be amplihed is also an issue, assumed in most cases as a 
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boundary condition. However, the consideration that in its initial state the Universe was 
homogeneous and isotropic to a high degree, excludes the presence of large scale magnetic 
helds. The seed magnetic field is therefore considered to be generated by some ‘battery 
mechanism’, the most common being the ‘Biermann battery’ (Biermann 1950), i.e. a 
mechanism for generating large scale electric currents based on the thermoelectric effect. In 
this respect it is worth reminding the reader that any battery mechanism depends on the 
possibility of imparting and sustaining different velocities to the electrons and protons in 
a plasma, thus creating an electric current. Obviously, in order for such an electric current 
to give rise to a steady (rather than transient) magnetic field it is necessary that the line 
integral of the proton-electron velocity difference along a closed, large scale circuit does not 
vanish. This leads directly to the notion that the generation of magnetic helds should be 
naturally associated with motions in which the integral /v-dl does not vanish, i.e. motions 
with non-zero circulation, a quantity also thought to be zero at the early Universe but 
generated later, in situations which involve the presence of dissipation, thus suggesting 
dissipation as a necessary condition for the generation of large scale magnetic helds too. 

Herein we revisit the theory of battery generated cosmic magnetic helds by employing 
a specihc mechanism for producing such currents in a specihc astrophysical setting, namely 
the Poynting-Robertson drag ehect on the electrons of an Advection-Dominated Accretion 
Flow (Narayan & Yi 1994, 95; hereafter ADAF) in the vicinity of an accreting object. We 
examine the time evolution of these magnetic helds in a variety of several astrophysical 
settings. We hnd that this ehect may be signihcant enough, at some phase in the evolution 
of the accreting object, as to not require any further amplihcation in order to account for 
the observed, present day values associated with the magnetic helds of the class of objects 
in question. 

In § 2 we give a brief outline of our mechanism and examine qualitatively its 
implications and ehects in diherent astrophysical settings. In § 3 we integrate numerically 
the associated equations and compare our results to those derived earlier on qualitative 
basis. Finally, in § 4 the results are summarized and conclusions are drawn. 


2. The model 

The model is simple: Consider a point-like astrophysical energy source of luminosity L, 
and an optically thin, fully ionized ADAF with rotational velocity = Arhlx at distances 
En ^ r < Tout around the above energy source (Dk is the Keplerian angular velocity, and 
A < 1 is a constant parameter of the how). Assuming that its scattering depth is sufficiently 
small that it scatters photons from the center only once, each plasma electron feels an 
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azimuthal force equal to 


Fp-R — 




( 1 ) 


A'lrr'^c c 

against its direction of rotation in the region between r[^ and rout where the plasma is 
optically thin. This is the Poynting-Robertson radiation drag. The positive ions feel a force 
{rrii/meY times smaller which can be neglected. As a consequence of the above force, the 
electrons tend to lag behind the protons, and an azimuthal electric current density 
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will develop^. Now, during the ADAF phase, the plasma is turbulent with a turbulent 
plasma conductivity at (see below), whereas when the ADAF phase ceases, the plasma 
conductivity attains a much higher value, uspitzer = 10^^ s“^ (e.g. Zombeck 1992; T„ 

is the plasma temperature in units of 10"’ K). On average, the effect of the hnite plasma 
conductivity cXeff can be parametrized as an extra azimuthal drag force equal to 


■^drag 
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( 3 ) 


acting on the electrons (similarly for the ions). The azimuthal component of the microscopic 
momentum (or force) equation for the electrons can be written as 


~ &{E --^P-R + -^drag — 0 , 


( 4 ) 


where E is the azimuthal component of the electric held; and Vr is the radial how velocity; e 
is taken positive. Combining eqs. 0.0. o. and (1) the axial component of the induction 
equation dB^/dt = —{c/r)d{rE)/dr yields 


dB^ 1 5 J La'pATtR 

dt r dr \ Ance 


rVrB;^ + 


dB., 


dTrcTeff dr 


( 5 ) 


^In order to keep the analysis simple, we will assume a cylindrical distribution for the magnetic field 
inside the geometrically thick ADAF, i.e. we will neglect its radial component and the associated magnetic 
field tension. 

"^The reader should not worry that ‘flux-freezing’ is not satisfied as electrons and ions move at slightly 
different azimuthal velocities. Electric currents generate magnetic forces which balance the forces Fg and 
Fi acting on electrons and ions respectively. The reader can convince him/herself that when the problem is 
solved, the two plasma constituents perform Fe,i x drifts across the above force fields, i.e. they 

do not move together even under conditions of flux-freezing (this analysis remains outside the main scope 
of the present work, and can be found in Ciolek & Mouschovias 1993). 
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in the range of radii (rin,rout)- ^ is defined through the relation /dr, and 27r\E' 

is the growth rate of the total magnetic flux contained within radius r. One should note 
that eq. (j^) includes both the advection of magnetic held through the term —WrE^, as well 
as its diffusion through the term jTna^isirdBjdr\ The connection between magnetic 
held and circulation discussed in the introduction, then, suggests that in an accretion disk, 
the turbulent magnetic dihusivity should also be proportional to the quantity ehecting the 
dissipation in the accretion disk, namely the turbulent viscosity parameter. Thus one may 
write (e.g. Reyes-Ruiz & Stepinski 1996) 


= cV(4’rCT.ff) ~ f’m!' ~ , 


( 6 ) 


where, a < 1 (the usual Shakura-Sunyaev (1973) a parameter) and < 1 (the ratio of the 
local how temperature to the virial temperature) are taken to be constant how parameters 
(we remind the reader that A = fl/flK = (1 — 5^/2 — a^^^/2)^/^; Narayan & Yi 1994, 95); 
Vm is the ratio of magnetic dihusivity to viscosity, the so-called magnetic Prandtl number. 
Note in addition that in an ADAF, the radial velocity can be related to the azimuthal one 
through the relation 

Vr ~ -a^rQx ~ -Ves/iVmr) . (7) 

Substituting all the above in the induction equation leads to 


OB, 

dt 
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r dr 
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The r.h.s. of eq. (|) consists of a source term (independent of defr), an advection term, 
and a dihusion term. The source term is the Poynting-Robertson azimuthal electric current 
and leads to held growth which is fastest at r = r^^. This gives rise to a distribution of 
closed poloidal magnetic held loops, with their center somewhere inside r^. As we will see 
in the next subsection, when VmA2, the magnetic flux contained inside r^^ grows almost 
linearly in time, as long as the ADAF persists. When rotation and accretion cease, and 
(Jeff ~ (Jspitzer, this flux wUl remain effectively trapped in the plasma over timescales 

2 

r ~ ——— ~ Ax'^MqTq^ yr . (9) 

^Spitzer 

Mn is the mass of the central object in units of 10” Mq-, x = r/r^ is the radial coordinate 
in units of the central object’s Schwarzschild radius rs. 

Obviously, this mechanism will cease being ehective when the magnetic held becomes 
strong enough to ahect the dynamics of accretion. In addition, a magnetic held which 
becomes dynamically important will also ahect the torques from the ADAF onto the star 
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(or vice versa). We prefer to defer discussion of the implications of this aspect of the 
problem and the rotational evolution of the compact object to a future publication; as 
indicated by the present analysis these could indeed be signihcant. 

It is of interest to examine the time scales necessary to accumulate enough flux via 
the above mechanism in the most optimistic case with Vm > 2, for the corresponding 
magnetic held to affect signihcantly the dynamics of accretion, as this may be an important 
component to the formation of jets, which are ubiquitous in accreting compact objects. 
Assuming an almost constant growth rate of magnetic hux d' of the order of the hrst term 
in the r.h.s. of eq. (^, 

^ Gcm\ (10) 

yyrs / 

where, / is the fraction of the source luminosity compared to the Eddington luminosity. One 
should note that the latter expression is independent of the mass of the system and depends 
only on its luminosity through the parameter / (for compact objects x is considered to be 
of order unity). However, the magnetic held itself is bounded by considerations associated 
with equipartition of the magnetic energy density with that of the ram pressure of almost 
free-falling matter, 

Heg ~ 9 • 10^ G , (11) 

where m is the accretion rate in units of the Eddington rate. Folowing Narayan & Yi 1994, 
we took a = 0.3, = 0.38, and A = 0.2. Assuming further that B{r < nn) ~ 'h(rin)/rj^, 

and that / ~ m/x, one obtains 

Teg ~ 4 • 10-V-2®Mo^-^m-°-5 yr . (12) 

This value is comparable to the Eddington time (~ 10® yrs) for an AGN with a massive 
central black hole (Mg ~ 1), but decreases very quickly with the decrease of the mass. 


3. Numerical Simulation 

The detailed time evolution of the magnetic held can be computed by direct integration 
of the time evolution equation (|^). We introduce the dimensionless variables 

r = r/rin , i = t/R , b = B/Bo , ip = 'iRo/{rf^Bo) , (13) 

where, to = {flK{rin)c(0~^'^^o = LaTA/(ATTcerfj^a^). The induction equation can now be 
rewritten in dimensionless form applicable in various astrophysical situations. 
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The induction equation is integrated numerically in the interval 1 < f < fout- In 
order to keep the number of free parameters to a minimum, we take fout 3> 1. We follow 
the numerical procedure described in the Appendix of Mouschovias & Morton (1991), and 
implement the ‘free boundary’ condition db/df = 0 at f = fin, allowing thus any amount of 
magnetic flux that would be advected inside f^ to enter unimpeded. We discovered that, 
for magnetic Prandtl numbers Vm > 2, the dimensionless magnetic flux inside dimensionless 
radius f = 1 rises almost linearly with time, with a dimensionless flux growth rate '0 of 
order unity. In fact, when Vm > 2, eq. (0) has the steady-state solution 


6(f) 


{2Vm - 1) 


4 ■ 

fO.5 ’ 


( 16 ) 


with, ip = 2{Vm — 2)/{2Vm — !)• We understand this interesting behavior as follows: The 
Poynting-Robertson effect generates an azimuthal electric current resulting in closed loops 
of poloidal magnetic held with a distribution of positive magnetic hux inside (where 
the magnetic held points, say, upward) and negative hux outside, i.e. in the regions with 
opposite held polarity. The rate of positive hux growth inside r-^ is given by the hrst term 
in the r.h.s. of eq. ( 0 - When Vm 'T. 2, dihnsion is dominant, and the distribution of 
negative hux outside rin dihuses outward to larger distances. As a result, positive hux is 
generated without limit (except of course for the equipartition limit) in the interior. On 
the other hand, when Vm < 2, the inward advection of the negative hux dominates over its 
outward dihnsion. As a result, the steady-state positive hux accumulated inside remains 
hnite, since its rate of growth by the Poynting-Robertson mechanism is balanced by the 
rate of inward advection of the return (negative) hux. 


In hgure 1 we exhibit the spatio-temporal evolution of the magnetic held and the 
associated hux resulting from our numerical solution to eq. ([0) . In hgure la we show the 
evolution of the radial distribution of the magnetic held as a function of the dimensionless 
radius r/vm for Vm = 2.5. Each curve corresponds to an increment in time by lOto- In 
hgure lb we exhibit the evolution of the magnetic hux held within as a function of time 
for diherent values of Vm- The reader should note that according to eq. (lA), for Prandtl 
numbers Vm > 2, a steady state situation, though it corresponds to a constant value for the 
magnetic held in the region r/vm > 1, demands a linear increase for the magnetic hux T for 
r/r-m < 1, with the return hux dihusing steadily to inhnity. 


In table 1, we introduce dimensions, and apply our mechanism to specihc astrophysical 
situations. The mechanism is most ehective at short scales, especially around compact 
objects (neutron stars, black holes of both stellar and galactic sizes). One could also suggest 
that magnetic hux generated quickly around active galactic nuclei, could be carried by 
winds out to galactic (and even extragalactic) scales (similar mechanisms have already been 
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proposed in order to account for the seed galactic dynamo field, e.g. Chakrabarti, Rosner, 
Vainshtein 1994). 


When we apply our mechanism to the most interesting stellar case, namely that of 
the Sun, we consider a scenario in which the innermost, optically thick, hydrogen burning 
core grows toward its main sequence values ~ O.lrQ, and ~ O.IMq respectively, or 
a; ~ 2 X 10®. We then consider that the remainder of the stellar mass (~ 1 Mq) accretes 
onto this component in the form of an ADAF, releasing and radiating away its gravitational 
energy. In this respect, the recent discovery of X-ray emission from objects in this class 
(Feigelson, private communication) is, to say the least, encouraging. Finally, we assume 
that the newly forming core radiates at some fraction / = 0.1 of the Eddington luminosity 
that corresponds to a mass equal to O.IMq (or equivalently L = 332 Lq; Casanova et 
al. 1995; Gross et al. 1998). The above are supposed to be crude approximations during 
the hnal stages of almost free-fall protostellar collapse. In the general protostellar collapse 
phase, 

T ~ 0.4Tq /Mif (rc/0.1 ro)-^-®(f/10® yr) . (16) 


Note that because of the linearity between the mass and the radius of the initial 
(proto)stellar cores, the above estimate, in accordance with eq. (pRD , is independent of 
their precise values, depending only on their ratio x. As one can see, for / ~ 0.1, it 
takes several million years to build the present day total observed solar dipolar magnetic 
flux Tq = 1 G • Tq. Considering the simplicity of these arguments, the agreement with 
observation is not unreasonable. 


We would like to discuss what is the geometry of the field thus generated and its fate 
after the mechanism ceases to be important (i.e. after the optically thin accretion flow 
dies out): The held is contained at the very core, in a layer around r^. In the absence of 
turbulence, the diffusion time-scale is of the order of several billion years, and the held will 



Un 

L 

M 

Uq 

Beq 

Black Hole 

3 • 10® 

10^ 

10 

7 ■ 10-^ 

5 • 

10^ 

Neutron Star 

10® 

10^ 

1 

10“® 

2 • 

10^ 

White Dwarf 

10® 

10 

1 

2 • 10^ 

2 • 

10^ 

Protosun 

7-109 

332 

10-1 

2 • 109 

8 • 

10® 

AGN 

3 • lO^® 

10^2 

10® 

7-10^ 

5 • 

10® 

Galaxy 

1022 

lO'o 

lOlo 

6 • 10®® 

6 • 

10-® 


Table 1: Tin in cm; L in Lq; M in Mq; Teq in yr; B^q = Teq/r?n G. 
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remain effectively trapped inside the solar radiation zone, i.e. below the turbulent surface 
convection zone. Any small fraction of this held which diffuses out of the sun will assume 
an overall dipolar geometry at large scales. In our mechanism the held is generated in situ, 
rather than being brought in from large distances. This magnetic held geometry should 
naturally interact with stellar diherential rotation, leading thus to torsional oscillations 
which might account for the solar cycle (e.g. Layzer, Rosner & Doyle 1979; Goode Sz 
Dziembowski 1991). 


4. Discussion 

The scenario described in this paper is not new. Previous researchers considered the 
rotation of solar coronal layers (Cattani & Sacchi 1966) as the source of the solar magnetic 
held, and the rotation of protogalaxies with respect to the cosmic microwave background 
radiation (Harrison 1970; Mishustin & Ruzmaikin 1972), as a means for producing the 
observed galactic held and (correctly) obtained very weak magnetic held values in both 
cases. However, to the best of our knowledge, this is the hrst time that this scenario has 
been applied to other classes of objects, in particular to optically thin accretion hows onto 
compact objects, where it appears to have its greatest impact, as the hgures in table 1 
indicate. 

Concerning the galactic magnetic helds, or better the value of the galactic magnetic 
hux Tg ~ 3 X 10^® G cm^ (this is mainly azimuthal and conhned on the plane of the galaxy) 
which presents the largest discord between the predictions of our model and observations 
(see table 1), one could conceivably consider its generation during an early AGN phase of 
our galaxy, instead of the era of galaxy formation; as such one could argue for the generation 
of the requisite hux near a massive black hole at the galactic center and its dihusion to 
large distances as discussed in the previous section, or as invoked by Chakrabarti, Rosner 
& Vainshtein (1994). Assuming that such an AGN stage lasted ~ 10® yrs, i.e. a few times 
the Eddington time (a quantity which is independent of the mass of the compact object), 
during which the mass galactic center black hole presumably grew to its present size and 
the associated magnetic held reached its equipartition value, eq.(0) yields a poloidal hux 
roughly 10^ times smaller than that observed. However, the galactic diherential rotation in 
the ensuing 10^® years (corresponding to roughly 100 galactic rotations at the outer regions 
of the galaxy and 100 - 300 times this value for its inner regions) could easily amplify this 
hgure to its presently observed value! 

Gonsidering in particular the generation of the helds in stars, this scenario can yield 
observationally interesting values for the magnetic held only if one is willing to consider 
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that during the pre-main sequence collapse of the cloud out of which the star formed, an 
advection-dominated accretion flow was probably present, and only if one is willing to 
accept values for the luminosity of this object (as given by /) and the duration of this 
process similar to those given above. However, even for lower values for / or for times 
shorter than ~ 10^ years, as in the above galactic case, the magnetic flux in the stellar 
interior could subsequently be amplihed linearly with time by the interior differential 
rotation to observationally interesting values. Given the very general scope of this work we 
would not like to pursue these issues beyond this level at present. We feel however that the 
assumptions about the pre-main sequence protostellar evolution employed in our arguments 
are not unreasonable. Interestingly enough, were the stellar magnetic held to be generated 
as discussed in the previous section, it should be strong and largely conhned to the deep 
stellar interior, leading naturally to the torsional oscillations between differentially rotating 
stellar layers which have been purported to account for stellar cycles. 

One cannot fail to notice that, when the outside magnetosphere can support electric 
currents (i.e. when there is no vacuum outside), the magnetic held geometry obtained 
in our numerical simulations, leads naturally to the electromagnetic interaction between 
the central region and the extended disk. The reason is that the diherential rotation 
‘twists’ the newly formed held, generating thus magnetic torques on the central part. The 
magnetic held twisting leads to another very important ehect, namely the overall opening 
of the closed loop geometry (Newman, Newman & Lovelace 1992; Goodson, Winglee & 
Bohm 1997). 1^ As one can see in the above references, the axial part of the expanding 
plasma + magnetic held remains well collimated. Under certain physical circumstances 
(Shibata, Tajima & Matsumoto 1992; Matsumoto, Matsuzaki, Tajima & Shibata 1997), 
this ehect might also proceed explosively, leading thus to the expulsion of a signihcant 
fraction of the disk material in the form of collimated fast axially moving plasma outhows 
(the astrophysical plasma gun, Gontopoulos 1995; Matsumoto et al. 1996). 

In this respect it is of interest that the time necessary for the magnetic held to become 
dynamically important is larger than the Eddington time for the most massive AGN with 
M>3 ■ 10® Mq. Thus one could speculate that the bright, optically selected QSOs have 
not had yet time to generate dynamically signihcant amounts of magnetic hux. One could 
further speculate that when they do so, the associated magnetic held would become an 
important component of the dynamics of accretion, facilitating the production of radio jets 
and leading to conversion of these objects from radio quiet to radio loud AGN. 


^In a simple ‘mechanistic’ interpretation, this can be thought of as the dynamical spring-like release 
through unwinding, of the azimuthal magnetic field generated by the differential rotation of the magnetic 
field footpoints. 
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The evolution of AGN and in particular its possible relation to their radio properties 
are far from being understood. The arguments put forward herein, associated with the 
generation and evolution of their magnetic helds may offer a new way of looking into this 
problem. Interestingly, the models of Falcke & Biermann (1995) which provide scalings of 
the radio loud objects from the quasar to the galactic binary source domain, do depend, 
to a large extent, on equipartition arguments such as those discussed above. It would 
be worthwhile to examine whether these arguments could be put in the context of an 
evolutionary scenario of both the mass or luminosity and the associated magnetic held. 

To conclude, it seems that the mechanism discussed above enters as a new, previously 
unaccounted for, player in the ongoing studies of pre-main sequence stellar collapse, the 
generation of stellar and galactic magnetic helds, the spin evolution of accreting magnetic 
stars, and the generation of magnetically driven jets. 
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signihcantly to its improvement. D.K. would like to acknowledge an informative discussion 
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Figure Captions 

Fig. 1. — Magnetic field evolution across the rotating plasma. In (a), we plot the 
magnetic field b = B^/Bg distribution as a function of the dimensionless radius r/rin, at time 
intervals lOtg, for Vm = 2.5. In (b), we plot the growth of the magnetic flux tjj = 
contained inside r = r\n with time, for different Prandtl numbers Vm- One sees clearly the 
fast growth of the held around r-m-, and the almost linear growth of the central magnetic 
hux for Vm > 2. The present (kinematic) integration should not be continued beyond Teg, 
when the held becomes dynamically important. 






